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ABSTRACT 

Aims. We investigate the environment in the vicinity of the supernova remnant (SNR) G59. 5+0.1 and identify all young stellar objects 
(YSOs) around the SNR, to derive the physical properties, obtain insight into the star-formation history, and further see whether SNR 
G59.5+0. 1 can trigger star formation in this region. 

Methods. We have performed the submillimeter/millimeter observations in CO lines toward the southeast of SNR G59.5+0. 1 with the 
KOSMA 3m Telescope. High integrated CO line intensity ratio R;co(3-2)/'co(i-i) identified as one good signature of SNR-MCs (molec- 
ular clouds) interacting system. To investigate the impact of SNR G59. 5+0.1 on the process of star formation, we used GLIMPSE 1 
Catalog to select YSOs (including class 1 and class II sources). 

Results. CO emission in an arc-like shape and mid-infrared 8.28 fim emission are well coincident with SNR G59.5+0. 1, which has 
the total mass of 1.1 x 10"* Mq and fully cover open cluster NGC 6823. Three molecular clumps were identified in the CO molecular 
arc, each clump shows the broad line wing emission, indicating that there are three outflows motion. The integrated CO line intensity 
ratio (R;coi,_ 

2i/'co(2-ii) f*-"^ whole molecular arc is between 0.48 and 1.57. The maximum value is 1.57, which is much higher than 
previous measurements of individual Galactic MCs. The CO molecular arc has a line intensity ratio gradient. SNR G59. 5+0.1 is in 
adiabatic expansion phase. The age of the SNR is 8.6 x lO"* yr. Based on GLIMPSE I Catalog 625 young stellar objects (YSOs) 
candidates (including 176 class I sources and 449 class II sources) are selected. The timescales for class 0, class I and class II sources 
are < lO'' yr, ~ 10' yr, and ~ 10* yr, respectively. The number of YSOs are significantly enhanced in the interacting regions, indicating 
the presence of some recently formed stars. 

Key words. ISM: clouds — ISM: individual (G59. 5+0.1) — ISM: molecules — stars: formation — supernova remnants 



] 1. Introduction 

. Massive star has significant impact on the morphology and 

■ chemical evolution of the surrounding Interstellar Medium 
[ (ISM) through UV-radiation, stellar winds, and supernova ex- 
. plosion. Stars form from the densest environments in molecu- 
' lar clouds (MCs), often clustered together spatially in groups 
] ranging from a few sources to many thousands. When a super- 
' nova explodes near MCs, shock generated by supernova remnant 
' (SNR) may compress some condensed MCs to collapse. On the 
, other hand, the shock can enhance abundances of the different 

■ molecular species with respect to quiescent cloud conditions. 
So the SNR-MCs interaction plays an essential part in the pro- 
cesses of star formation and the evolution of ISM (Reynoso et al. 
12001) . These associations are established by the detection of OH 
1720 MHz masers, the broad emission lines, high ratio of '^CO 
7 = 2-1 to 7 = 1-0 integrated line intensity, y-ray emission 
from MCs toward SNR as well as the morphological relation be- 
tween SNR and MCs, and so on (Frail et al. 119961 Green et al. 
[19971 Seta et al. 1998. Wilner et al. 19981 Hua ng & Thaddeus 
[19861 Byun et al. [20061 Hewitt & Yusef-Zadeh [^09l ). Jiang et 
al (|2010) summarized the criteria of the SNR-MCs association. 
In addition, some results from numerical simulation have specu- 
lated that star formation can be triggered by the SNR-MCs inter- 
action (Melioli et al. i2006i & Leao et al. l2009l l. However, only a 



few cases about the triggered star formation have been proposed 
from the observed results (Junkes et al. |1992t Parons et al. 120091 
Xu et al 12011 al l . 

SNR G59.5+ 0.1, a t a distance of 2.1-2.3 kpc (Xu et al. 
2005; Billot et al. 120 101 ). is situated in the direction of Vulpecula 
OB at the position {I,b) - (59.58, 0.12). According to Xu et 
al. (2005), the age of SNR G59.5+0.1 is ranges from 10^ to 
10"^ years old. Taylor et al. ([1992]l described this object as a 
shell-type SNR with nonthermal spectral index (a < -0.4) and 
with a diameter of 15'. They also mentioned that the SNR is 
possibly associated with the HII region Sh2-86. Billot et al. 
(2010) found a circular structure centered at the position of the 
SNR with a diameter of 15', which confirms the value found 
by Taylor et al. (1992) . Furthermore, Billot et al. (120 lOD noted 
a high density of YSOs lined up along two arcs at the north- 
ern and southern rims of SNR G59. 5+0.1, suggesting that the 
interaction of the expanding shell of SNR G59. 5+0.1 and the 
neutral gas of Sh2-86 could have triggered the formation of 
young stars in the cluster. NGC6823 is an open cluster located 
at R.A(J2000)=19M3'"08^4, Dec(J2000)=23°18'00" wit h a di- 
ameter of ~ 16' (Kharchenko et al. 2005). Riaz et al. (120121) 
suggested that NGC6823 have possibly experienced a recent 
burst of star formation, which may be caused by a supernova 
explosion of massive O star However, they suggested that SNR 
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Fig. 1. CO J - 2-1 intensity map (black contours) integrated 
from 1 8 to 35 km s"' , overlayed on the the mid-infrared 8.28-//m 
MSX emission map (color scale). The contour levels are 30, 40, 
... , 90% of the peak value (83.9 K km s"'). The red star repre- 
sents the center of SNR G59. 5-1-0.1 and its extent is outlined by 
a red dashed circle. The radius of SNR G59.5+0. lis 15' (Taylor 
et al. 1992 & Billot et al. 2010). The white box outlines the field 
seen on the Fig.2. 



G59. 5-1-0. 1 is not old enough to trigger any of their identified star 
formation. 

Motivated by the supposed association of SNR G59. 5+0.1 
with molecular gas and with possibly triggered star formation, 
we have performed '^CO J = 2-1, '^CO 7 = 3 - 2, and 
'^C07 = 3-2 observations toward G59.5+0.1. The observa- 
tions are described in §2, and the results are presented in §3. In 
§4, we discuss how our data lend supportive evidence of the trig- 
gered star formation in the interacting region. The conclusions 
are summarized in §5. 



2. Observations 

The mapping observations of G59. 5-1-0.1 were made in '^CO 
/ = 2 - 1, i^co 7 = 3 - 2 and '^CO 7 = 2-1 lines using 
the KOSMA 3m telescope at Gornergrat, Switzerland, in April 
2004. The half-power beam widths of the telescope at the ob- 
serving frequencies of 230.538 GHz, 345.789 GHz and 220.399 
GHz are 130", 80", and 130", respectively. The pointing and 
tracking accuracy is better than 20". The accuracy of the ab- 
solute intensity calibration to be better than 15% (Sun et al. 
12008 ). The DSB receiver noise temperature was about 170 K. 
The medium and variable resolution acousto optical spectrome- 
ters have 1501 and 1601 channels, with total bandwidth of 248 
MHz and 544 MHz. The channel widths of 165 and 340 kHz 
correspond to velocity resolutions of 0.21 and 0.29 km s"', re- 
spectively. The beam efficiency Beff is 0.68 at 230 GHz and 220 
GHz. The beam efficiency Beff is 0.72 at 345 GHz. The for- 
ward efficiency F^ff is 0.93. Mapping observations are centered 
at RA(J2000)=19H3™47.85^ DEC(J2000)=23°21 '26.34" using 
the On-The-Fly mode, the total mapping area is 24' x24' in '^CO 
7 = 2 - 1 and '^CO 7 = 3-2 with a 1' x 1' grid, while the map- 
ping area in '^CO 7 = 2 - 1 is 10' x 15'. 
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Fig.2. '^CO 7 = 2-1 intensity map (gray scale), overlayed on 
'^CO 7 = 3-2 (blue contours) and CO 7 = 2 - 1 (white con- 
tours) intensity maps. The contour levels of each CO molecule 
are 30, 40,..., 90% of the peak value. The '^CO 7 = 3 - 2 and 
'^CO 7 = 2-1 peak values are 79.0 and 15.5 K km s ', respec- 
tively. Letter A, B and C indicate the different MC clumps. The 
dot symbols mark the mapping points of '^CO 7 = 2-1 and 
'^CO 7 = 3-2, while the black rectangle outlines the mapped 
field of '^C0 7 = 2- 1. 

The data were reduced using the GILDAS/CLASS Q pack- 
age. The correction for the line intensities to the main beam tem- 
perature scale was made using the formula Tmb = (^'eff/Beff x 

3. Results 

3.1. Molecular emission 

Figure 1 presents the emission map (color scale) at 8.28 jum. 
The emission shows a half-shell morphology, which is associ- 
ated with SNR G59.5+0. 1 with a diameter of 15'. The integrated 
intensity map of '^CO 7 = 2 - 1 is overlapped the 8.28 fim emis- 
sion located in the white boxed area of Figure 1, which cov- 
ers open cluster NGC6823 marked in a green dashed circle (see 
Figure 8). The '^CO 7 = 2-1 emission coincides well with the 
8.28 fxm emission at the southeast. In order to further analyze 
the morphology of CO molecular gas, we make the integrated 
intensity map of '^CO 7 = 3 - 2 and '^CO 7 = 2-1 (Figure 
2). In Figure 2, '^CO 7 = 2-1 emission is coincident with that 
from I^CO 7 = 3 - 2 and '^CO 7 = 2 - 1. The '^CO 7 = 2-1 
and '^CO 7 = 3-2 emission appear to show a arc-like structure, 
while the '^CO 7 = 2-1 emission may trace the compact core 
of molecular gas. We find three cloud clumps in this structure, 
each clump is designated alphabetically, clump A, clump B and 
clump C. Because the mapping area in '-^CO 7 = 2 - 1 is smaller 
(see section 2), we did not detected the '-^CO 7 = 2-1 emission 
from the clump C. Figure 3 shows '^CO 7 = 2-1, '^CO 7 = 3-2 
and '^CO 7 = 2-1 spectra averaged over clump A, clump B and 
clump C, respectively. Line profiles of '^CO 7 = 2-1 and '^CO 
7 = 3-2 for each clump appear to be broadened, the velocity 
components are mainly located in interval 18 to 35 km s The 

' http://www.iram.fr/IRAMFR/GILDAS/ 
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Fig. 3. Averaged spectra of molecular line '^CO 7 = 2-1, '^CO 7 = 3-2 and '^CO 7 = 2-1 over each clump. The dashed lines 
mark the position at the zero intensity. 



Table 1. Observed parameters of each line 



Name 




'-CO 2-1 






'^CO 3-2 






'^CO 2-1 






Tmb 


FWHM 


Vlsr 




FWHM 


Vlsr 




FWHM 


Vlsr 




(K) 


(kms-') 


(km s-') 


(K) 


(km s-') 


(km s-') 


(K) 


(kms-') 


(km s-') 


Core A 


9.9 (0.6) 


3.7 (0.2) 


27.0 (0.1) 


14.5 (0.4) 


3.9(0.1) 


27.0 (0.1) 


5.4 (0.2) 


2.7 (0.1) 


27.3 (0.1) 


Core B 


9.2 (0.2) 


4.4 (0.1) 


27.4 (0.1) 


9.4 (0.2) 


4.4(0.1) 


27.6 (0.1) 


3.6 (0.2) 


3.6 (0.2) 


27.6 (0.2) 


CoreC 


13.9 (0.3) 


4.4 (0.1) 


28.8(0.1) 


18.3 (0.3) 


4.0(0.1) 


29.0 (0.1) 









Note: The position of clump A is at RA(J2000)=19''43"'51.85'^ and DEC(J2000)=23°28'26.34"; The position of clump B 
is at RA(J2000)=19M3'"47.85^ and DEC(J2000)=23° 19'26.34"; The position of clump C is at RA(J2000)=19M3"15.85^ and 
DEC(J2000)=23°13'26.34". 



spectral profile of clump A may show a broad outflow wing and 
has two weaker components at 23 km s"' and 31 km s"'. We 
have made Gaussian fits to all the spectra. The fitted results are 
summarized in Table 1. From Table 1, we derive systemic veloc- 
ities of ~27.1 km s"', ~27.5 km s"\ and ~28.9 km s"' in clump 
A, clump B, and clump C, respectively. According to the galactic 
rotation model of Fich et al. (1 19891 1 together with Rq = 8.5 kpc 
and Vq - 220 km s"' , where V© is the circular rotation speed of 
the Galaxy, we obtain a kinematic of 2.3 kpc for the distance 
of all clumps, which is roughly consistent with the photometric 
distance (2.1 ±0.1 kpc) of NGC 6823 (Guetter 1992; ). 

Assuming local thermodynamical equilibrium (LTE) and us- 
ing the '^CO 7 = 2-1 line, the column density of the clumps in 
cm - is estimated as (Garden et al. ll991l Xu et al.|2010) 



clu 



1.08 X 10' 



, (Fex + 0.92) 
exp(-16.62/rex) 



rrnbX 



rdv 



[1 -exp(-T)] 



.(1) 



where dv is the velocity range in km s"', Tex is the excita- 
tion temperature in K, and t is the optical depth for '^CO 
7 = 2-1 line. Tex is estimated following the equation Tex = 



ll.l/ln[l + l/(rmb/ll-l + 0.02)], while we calculate t accord- 
ing to 



r^bC'CO) 1 - exp(-T) 



r„b(i3CO) l-exp(-T/89)' 



(2) 



where T^b is the corrected main beam temperature. Here we 
assume the solar abundance ratio of ['^CO]/['-'CO]=89 (Lang 
1980; Anders & Grevesse 1989; Garden et al. fT99TT l. and that 
the '^CO 7 = 2-1 emission is optically thick. Additionally, 
we use the relation A^h, ~ lO^^A^nco (Dickman [19781 ). If the 
clumps are approximately spherical in shape, the mean number 
density H2 is n(H2)=1.62 x lO^^'^Nui/L, where L is the clump 
diameter in parsecs (pc). Furthermore, their mass is given by 
MHj=n(H2) ^ttL Vg'«(H2) (Garden et al. 0.991), where jt/^=1.36 
is the mean atomic weight of the gas, and ;«(H2) is the mass of 
a hydrogen molecule. The derived column density, mean num- 
ber density, and mass of each clump may be in error by as 
much as a factor of three depending on the exact value of the 
['2CO]/['^CO] abundance ratio (Garden et al. [19911 ), which are 
all listed in Table 2. 
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Fig. 4. P-V diagram constructed from CO 7 = 2-1 transition for clump A-C. (a)-(c) panels: Contour levels are 10, 20,..., 90% of 
the peak value, with a cut through the peak positions of each clump along the north-south direction, the peak positions are shown 
in the upper right corner, (d)-(f) panels: Contour levels are the same as above panels, with a cut along the east-west direction. The 
blue and red vertical dashed lines mark the velocity ranges of the blueshifted and redshifted emission for each clump, respectively. 
The offsets given in each panel are relative to RA(J2000)=19H3'M7.85^ and DEC(J2000)=23°21'26.34". 



Table 2. The physical parameters of the clumps in LTE. 



Name 


T 

' ex 


L 






M 




K 


pc 


(cm-2) 


(cm-3) 


(10^' Mo) 


Clump A 


15.0 


5.4 


1.4x10^^ 


4.2x10^ 


6.3 


Clump B 


14.3 


4.7 


0.9x10^2 


3.3x10^ 


3.2 


Clump C 


19.2 


6.3 


2.2x10" 


5.6x10^ 


1.4 



To determine the velocity components and confirm the out- 
flow, we have made position-velocity (PV) diagrams with the 
cuts through the peak positions of each clump along the north- 
south and east-west directions, as shown in Figure 4. From 
Figure 4, we can clearly identify three clumps. Each clump 
shows the bipolar structure marked by the blue and red vertical 
dashed lines. The blueshifted and redshifted components have 
obvious velocity gradients, in particular for clump A and clump 
B presented in panel (d) and (e) of Figure 4. The distributions of 
redshifted and blueshifted velocity components are a indication 
of outflow motion. In panel (d), we also identify a velocity com- 
ponent from 21 to 24 km s"', which is not from the component 
of the outflow, but related to a '^CO J - 2-1 weaker component 
peaked at -23 km s"' as seen from the '^CO J - 2 - I spectra 
of clump A in Figure 3. Using the velocity ranges obtained from 
the PV diagram, we made the velocity integrated intensity maps 
superimposed on each clump emission map, respectively. The 
distribution of redshifted and blueshifted velocity components 
in Figure 5 provides us further evidence for the bipolar outflow 
in each clump. 



o 23°24'00" 



Q 23°16'00" 
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Fig. 5. The velocity integrated intensity maps of '^CO J - 2-1 
outflows (red and blue contours) overlaid with the '^CO J - 2-1 
emission of each clump (gray scale). The red and blue contour 
levels are 20, 30,. ..,100% of the peak value. The red dashed arc 
represents SNR G59.5+0.1. 



When a supernova expands into the surrounding molecular 
clouds, the shocks resulting from the interaction of SNR with 
MCs can heat the surrounding gases. As the temperature of gases 
increases, the line opacities decrease as the upper J levels be- 
come more populated, then the integrated intensity ratio of '^CO 
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Fig. 6. '^CO 7 = 2-1 intensity map are superimposed on the 
Une intensity ratio map (color scale), the line intensity ratios 
C^icoo-D/icor-t)^ range from 0.48 to 1.57 by 0.16. The wedge in- 
dicates the line intensity ratios scale. 



7 = 3 - 2 to '^CO / = 2 - 1 (R, 



CO(3-2)/'COa-l) 



) can indicate the 



shock (Xu et al. l201 lall201 Ibk In order to obtain the integrated 
intensity ratio of '^CO 7 = 3 - 2 to '^CO 7 = 2 - 1, the 80" 
resolution of '^CO 7 = 3-2 data is convolved with an effective 

102". We calculated the integrated 
2 line in the same velocity range as 



beam size of Vl302-802 
intensities for '^CO 7 = 3- 

for '^CO 7 = 2-1 line. The integrated range is from 16 to 36 
km s"'. Figure 6 shows the distribution of R/co(3-i)//co(2-i) (color 
scale) overlaid with the distribution of '^CO 7 = 2-1 line inte- 
grated intensity (contours). The red dashed arc represents SNR 
G59.5-(-0.1. In Figure. 6, the whole arc-like molecular gas has a 
ratio value gradient increasing from northeast to southwest, sug- 
gesting that SNR shock is expanding into clump A-C and start to 
compress each clump. The ratio values for the clump A, clump 
B, and clump C are between 0.48 and 1.57. The maximum value 
is 1.57, which is much higher than typical value (0.55) for MCs 
in the Galactic disk (Sanders et al. 1 19931 1 and value (0.69) for the 
normal MCs in M33 (0.69 Wilson et al. ll997l l. and even higher 
than value (0.8) in the starburst galaxies M82 (Guesten et al. 
fT993T i. 
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Fig, 7, GLIMPSE color-color diagram [5.8]-[8.0] versus [3.6]- 
[4.5] for sources. The regions indicate the stellar evolutionary 
stage as defined by Allen et al. (2004). Class I sources are proto- 
star with circumstellar envelopes and class II are disk dominated 
objects. 



sources are protostars with circumstellar envelopes (polygon), 
class II sources are disk-dominated objects (rectangle), and other 
sources. Using this criteria, we find 176 class I sources and 449 
class II sources. Here class I and class II sources are chosen to 
be YSOs. 

Figure 8 shows the spatial distribution of both class I and 
class II sources. From Figure 8, we note that class I and class II 
sources are not symmetrically distributed in the whole selected 
region, and are mostly concentrated in the north and southeast 
around SNR G59.5H-0.1, which is similar to the results of Billot 
et al. ( 120101 ). Regarding the geometric distribution of class I and 
class II sources, we can plot the map of the star surface density. 
Because class I sources is younger than class II sources and the 
spatial distribution of class II sources is similar to that for class I 
sources, we only plot the map of the star surface density for class 
I sources. This map was obtained by counting all class I sources 
with a detection in the 3.6 //m, 4.5 //m, 5.8 //m, and 8.0 fim bands 
in squares of 4' x 4', as shown in Figure 9 . From Figure 9 we 
can see that there are clear signs of clustering toward the region 
where SNR G59. 5+0.1 is close approximately to the surround- 
ing ISM. The existence of class I sources may also indicate star 
formation activity. 



3.2. Infrared emission 

To further search for primary tracers of the star-formation ac- 
tivity in the surrounding of G59.5H-0.1, we used the GLIMPSE 
I Catalog which consists of point sources that are detected at 
least twice in one band. The GLIMPSE I survey observed the 
Galactic plane (65° < |/| < 10° for \b\ < 1°) with the four mid-IR 
bands (3.6, 4.5, 5.8, and 8.0 fim) of the Infrared Array Camera 
(IRAC; Fazio et al. 120041) on the Spitzer Space Telescope. 
From the database, we selected 29469 near-infrared sources 
within a circle of 24' in radius centered on R.A.=19H2"50^ 
(J2000), Dec=-H23°30'00" (J2000). The size of this region com- 
pletely covers the extension of G59.5H-0.1 and the open cluster 
NGC6823. Figure 7 shows the [5.8] - [8.0] versus [3.6] - [4.5] 
color-color (CC) diagram. The regions in the figure indicate the 
stellar evolutionary stages based on the criteria of Allen et al. 
(I2004l l. Parons et al (|fD09l l, and Petriella et al (I2010al) . The 
near-infrared sources are classified into three regions: class I 



4. Discussion 

4.1. SNR G59.5+0.1 and MCs 

Because of the high confusion level in the surrounding of Sh2- 
86 structured emission, both Reach et al. (2006) and Billot 
et al. ([20T0t did not detect SNR G59.5+0.1. Similarly, we 
could not detected the SNR on the 8.28 /urn image, but de- 
tected the 8.28 fj.m emission which present the half-shell struc- 
ture around the SNR. The half-shell emission may be produced 
by SNR G59. 5+0.1 and/or the stellar winds of its progenitor 
We have performed the submillimeter/millimeter observations 
in CO lines toward the southeast of SNR G59.5+0.1. The '^CO 
7 = 2-1 molecular emission is coincident very well with the 
8.28 fim emission, which covers the whole size of NGC 6823. 
In addition, the '^CO 7 = 2-1 molecular emission shows the 
arc-like morphology. We identified three clumps in the molec- 
ular arc. Three bipolar outflows are detected in these clumps 
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Fig. 8. Positions of Class I and II sources relative to the '^CO 
7 = 2-1 MCs (grey scale) around SNR G59.5+0.1. The Class 
I sources are labeled as the blue dots; The Class II sources are 
labeled as the red dots. The black and green dashed circles rep- 
resent SNR G59.5+0. 1 and open cluster NGC6823, respectively. 
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Fig. 9. Stellar-surface density map of Class I candidates are su- 
perimposed on '^CO 7 = 2-1 MCs (grey scale) around SNR 
G59.5-I-0.1. Contours range from 2 to 10 stars (4arcmin)"^ in 
steps of 4 stars (4arcmin)"^. Icr is 0.7 (4arcmin)"^ (background 
stars). 



around SNR G59. 5+0.1 for the first time. The line profile of 
'^CO 7 = 2 - 1 in clump A appear to be broadened, which 
may be not only caused by the outflows, but also by shock from 



SNR G59.5+0.1. The maximum value of R 



/cO(3-2)//cO(2-l) 



1.6 



in the molecular arc, which is higher than that in the previous 
measurement of individual Galactic MCs. High R/co(3-i)//co(^-i) 
is suggested as a signature of the SNR-MC interaction system 
(Jiang et al. |2010[ Xu et al. 2011a). Together, these observations 
strongly indicate that the MCs are interacting with G59.5+0.1. 
SNR G59.5+0.1 just start to interact with the surrounding MCs, 



so the MCs have the R/co(3-2)//co(2-i) gradient along the shock di- 
rection. We suggest that the distribution of R/co(3-2)//co(^-i) value 
is a good signature of the SNR-MCs interaction system. The 
integrated CO intensity ratios (R/co(3-2)//co(^-i)) be used to 
study other triggering mechanism of star formation, such as the 
shocks of cloud-cloud collision, HII regions and galactic den- 
sity waves. From Table 2, the total mass of MCs associated with 
SNRG59.5+0.1 is 1.1 x 10^ M^. 



4.2. Recent star formation in tine MCs 

To summarize the results from investigating the spatial distribu- 
tion of YSOs, we clearly see that the distribution of YSOs (class 
I and class II sources) is clustered around the border of SNR 
G59.5+0. 1 . The best correlation is found with the CO molecular 
arc in the southeast where most of YSOs may belong to open 
cluster NGC6823. It is unlikely that they all are just foreground 
and background stars. It is more likely that those YSOs are phys- 
ically associated with the interacting regions between G59.5-H0. 1 
and MCs. Riaz et al. (2012) suggested that NGC6823 have pos- 
sibly experienced a recent burst of star formation, which may be 
caused by a supernova explosion of massive O star Also, these 
YSOs distribution is clustered together and shows a shell-like 
structure around G59.5+0.1. SNR G59. 5-1-0.1 is close approxi- 
mately to NGC6823. We suggest that SNR G59.5+0.1 possibly 
have triggered these YSOs formation and star formation in this 
cluster Class I sources occur in a period on the order of ~ 10^ yr, 
while the age of class II sources is ~10^ yr (Andre & montmerle 
1994). 

In an inhomogeneous medium, an SNR can be undergoing 
different evolutionary stages in different places at the same time. 
If an SNR is in radiative expansion phase, the radius of an SNR is 
-30 pc (HeilesT964). Here the radius of SNR G59.5-H0.1 is -10 
pc, the SNR has not yet been achieved the radiative phase, but 
in adiabatic expansion. The age of G59. 5+0.1 given by Sedov 
equation (Reynoso & Mangum l2001l l in units of 10"^yr: 



(3) 



where £51 is the remnant energy, we adopt a canonical value of 
10^' ergs. Rs is the remnant radius, while no is the the ambient 
number density in cm^^, which is determined by Table 2. Then 
we obtain that the age of G59. 5+0.1 is 8.6x10"* yr. The age of 
G59.5+0.1 are shorter than that of these YSOs, which is 10^ to 
10^ yr Hence the formation of the YSOs may not be triggered 
by the shock of SNR G59.5+0.1. Because SNR G54.4-0.3, SNR 
G24.7+0.6 and SNR IC443 are not old eno ugh, J unkes et al. 
( ri992D , Petriella et al (2011), and Xu et al. (1201 lai) suggested 
that the YSOs around the remnants may not be triggered by the 
remnants, but triggered by the stellar winds of the remnants pro- 
genitors. SNR G59.5+0. 1 is located in the direction of Vulpecula 
OB at the position {l,b) = (59.58, 0.12), hence the spectral type 
of SNR G59. 5+0.1 progenitor may be O or B. In the course 
of high-mass stars evolution the stellar winds go through three 
phases (Dwarkadas l2007l) . which last for about 4.6x10^ yr Such 
periods of stellar winds are sufficient to form YSOs. We con- 
clude that the formation of the YSOs may be triggered by the 
stellar winds from the high-mass progenitor of G59.5+0. 1 . 

Outflow is a strong evidence of the earlier star forming ac- 
tivity. We detected three outflows around SNR G59.5+0.1. The 
dynamic timescale of each outflow is given by equation t = 9.78 
X lO^'R/V (yr), where R in pc is the outflow size defined by the 
average of the radius of the blueshifted and redshifted lobes, V 
in km s"' is the maximum flow velocity relative to the cloud 
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systemic velocity. The obtained average dynamical timescales 
of outflow A-C are -7.4x10^ yr, -5.3x10^ yr, and -5.2x10^ 
yr, suggesting that there are three Class protostars (< 10^ yr). 
These Class protostars are distributed around SNR G59.5+0.1. 
SNR G59.5+0.1 with a age of 8.6x10"* yr may trigger the for- 
mation of these Class protostars. The observations with higher 
spatial resolutions are needed to further analyze these protostars. 

5. Conclusions 

We have presented the ^^CO 7 = 2-1, ^^CO / = 3 - 2, and 
'"'CO 7 = 2-1 molecular and infrared observations towards 
SNR G59.5+0.1. These results can be summarized as follows: 

1. The arc-like morphological association with SNR 
G59.5+0.1, the broadened emission Hnes, and the high 
integrated CO Une intensity ratio (R/co(3-2)//co(2-i))(~ 1-58) 
suggest that these cloud clumps are interacting with 
G59.5+0.1. The age of the SNR is 8.6 x 10"* yr The 
whole molecular gas in the southeast have a line intensity 
ratio gradient along the direction of shock, implying that 
shocks have driven into cloud clumps. We suggest that high 
R/co(3-2)//co(2-i) is identified as one good signature to study 
the triggering mechanism of star formation, such as the 
shocks of cloud-cloud collision, HII regions, and Galactic 
density waves. The cloud clumps have the total mass of 
1.1 X W Mq. 

2. The selected young stellar objects (YSOs) (including class 
I and class II sources) are concentrated and grouped around 
the interacting regions. It provides strong signpost for ongo- 
ing star formation. Comparing the age of G59.5+0.1 and the 
timescales of the stellar winds of G59. 5-1-0.1 progenitor with 
the characteristic star-formation timescales, we conclude 
that the YSOs may not be triggered by SNR G59. 5+0.1, but 
triggered by the stellar winds of G59.5-H0.1 progenitor. 

3. Three young outflows are detected in the clumps around 
SNR G59.5+0.1. It may be the first SNR around which find 
such more outflows in all the Galaxy SNRs. Taking into ac- 
count the age of SNR G59. 5+0.1 and class source, we 
find that SNR G59.5+0.1 may trigger the formation of these 
Class sources. It also may provide us a direct evidence for 
star formation triggered by SNR. The results of the present 
work confirm the sequential star formation on the basis of 
the different ages of star formation associated with SNR 
G59.5+0.1. 
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